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ABSTRACT 

Upflows observed at the edges of active regions have been proposed as the source of the 
slow solar wind. In the particular case of Active Region (AR) 10942, where such an upflow 
has been already observed, we want to evaluate the part of this upflow that actually remains 
confined in the magnetic loops that connect AR 10942 to AR 10943. Both active regions were 
visible simultaneously on the solar disk and were observed by STEREO/SECCHI EUVI. Using 
Hinode/EIS spectra, we determine the Doppler shifts and densities in AR 10943 and AR 10942, 
in order to evaluate the mass flows. We also perform magnetic field extrapolations to assess 
the connectivity between AR 10942 and AR 10943. AR 10943 displays a persistent downflow in 
Fe XII. Magnetic extrapolations including both ARs show that this downflow can be connected to 
the upflow in AR 10942. We estimate that the mass flow received by AR 10943 areas connected to 
AR 10942 represents about 18% of the mass flow from AR 10942. We conclude that the upflows 
observed on the edge of active regions represent either large-scale loops with mass flowing along 
them (accounting for about one-fifth of the total mass flow in this example) or open magnetic 
field structures where the slow solar wind originates. 



Subject headings: methods: data analysis — Sun: atmospheric motions 
radiation — techniques: spectroscopic 



Sun: corona - Sun: UV 



1. Introduction 

The Sun interacts with the whole heliosphere, 
and in particular with the planets of the solar sys- 
tem, through the solar wind. As the plasma /3 
(ratio of the plasma pressure to the magnetic pres- 
sure) is low in the low corona, the dynamics of the 
plasma is dominated by the magnetic field (frozen- 
in condition) implying that the plasma material 
is flowing along magnetic field lines. In particu- 
lar, the fast solar wind follows open magnetic field 
lines from solar coronal holes to the interplanetary 



space but the sources of the slow wind remain an 
open issue. 

Fast and slow winds can be distinguished ac- 
cording to their speeds (around 600 and 300 km • 
s _1 respectively) and their composition, but this is 
not the only difference between both these types. 
The fast wind is ex pelled from coronal holes, espe- 
cially at the poles ( Krieger et al.lll973 ) , and possi- 
bly from the intersections of chro mospheric net- 



1 Univ Paris Sud, Institut d'Astrophysique Spatiale, 
UMR8617, 91405 Orsay, France 

2 CNRS, Institut d'Astrophysique Spatiale, UMR8617, 
91405 Orsay, France 

3 Jeremiah Horrocks Institute, University of Central 
Lancashire, Preston, Lancashire, PR1 2HE, UK 



work boundaries in coronal holes (lHassler et al 
Il999h . The slow solar wind is not as well under- 
stood as the fast one, for various reasons probably 
related to its time variability. Its transient na- 
ture and its relation with large-scale coronal struc- 
tures, have been r evealed from both in-situ (e.g. 
Kilpua et al. 2009[ ) and remote-sensing obs erva- 



tions (e.g. the blobs of Sheelev et al. 1997t ). As 
for the coronal sources, the edges of coronal holes 
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(Coronal Hole Boundaries) have been proposed as 
the location of reconnection between coronal hole 
(CH) and non-CH magnetic fields, because of the 
differential rotation between these two kinds of 



2005; Wang fc Sheelev 



regio ns (jSchwadron et al 
2004). Recent spectrocopic and imaging observa- 
tions with SUMER/SOHO and XRT/Hinode seem 
to support this mechanism (jMadiarska et al . 2004, 



Subramanian et al. 2010). However, other loca- 



tions and mechanisms ha ve been proposed, such 



as streamer boundaries (lAbbo et al.l 120101) and 



the edges of Active Regions ( Koiima et al.lll999t 
Liewer et all 12004k iKo et al.l l2006h . As "we shall 



see below, this latter possibility has been very re- 
cently put forward in the context of an Active Re- 
gion close to an "open field" region, an issue we 
focus on in this Paper. 

Active regions (ARs) in the Sun's atmosphere 
are composed of closed multi-temperature loops 
in the solar corona. Recently a specific flows 
distribution has been shown for some ARs. It 



corre s ponds to redshi f ts ins ide loops (|Del Zanna 
2008L iTripathi et"aH l2009h and blueshifts at 



the edge of active regions (Doschek et al. 2007L 



Sakao et al.l 120071 lHarra et al l 120081 IBaker et al 



2009, iDel Zanna et al.ll201ll) . In the last two pa- 



pers, according to magnetic field extrapolations, 
blueshifts are observed along fanning out, far- 
reaching or even open field lines. Then the ob- 
served flows seem to occur at the boundaries be- 
tween active regions and coronal holes. Thus the 
outflows can supply mass to the solar wind, as 
suggested by scintillation measurements at 2.5 
solar radii, these outflows coming from an actu- 
ally open region at the edge of the active region 
(jKoiima et al.lll999h . 

Active Region 10942 has already been exten- 
sively studied: Fast upflows have been observed at 
the North-East edge of this AR from Hinode/EIS 



Dopp l er shifts in Fe XII 195.12 A (lHarra et al 



2008; IBaker et all l2009t iMcIntosh fc De Pontieu 



2009 ) and from a pparent flows in Hinode/XRT 



AR 10942 was connected by large scale loops to 
a magnetic dipole, approximately 400 arcsec away, 
which actually is AR 10943. This is a clear ev- 
idence for magnetic connection, but matter ex- 
change has never been quantified. Large scale 
loops connecting two active regions have been 



observed since Sky lab (Vaianal 



1976) and more 



recently with SDO (|Schriiver fc Title! 1201 if ). The 
connection between far-distance active regions has 
also been observed in the case of transequatorial 
l oops whether they are elo ngated stable structures 
( Farnik et al.l 1 1999L 12001 ) or epheme ral loops re 



lated to flare and filament eruption (jWang et al 



2007). Of special interest here are the spectro- 



scopic observations of a transequatorial loop which 
indicate that the loop plasma was multithermal 
and cover ed rough l y 2 or ders of magnitude in tem- 
perature (jBrosiusI 120061 ) . Moreover line-of-sight 



steady flows of the order of 30 to 40km.s~ 1 were 
detected and interpreted as a necessary condition 
for maintaining the loop structure. 

The above discussions on flows in active re- 
gions should not lead us to forget the issue of the 
"rest" wavelengths used to define absolutely the 
flows, even if these flows are relatively important 
in AR. Because the region taken as a reference 
is often the quiet Sun (QS), the issue of average 
temperature-d e pende nt flows in the QS is critical. 



Dadashi et al.l (|20111) mention average line shifts 
at 1 MK < T < 1.8 MK bluer than those ob- 
served at 1 MK (about -1.8 ± 0.6 km.s -1 ) trans- 
lating into a maximum Doppler shift of -4.4 ±2.2 
km.s -1 around 1.8 MK. If we assume that the ac- 



one 



dSakao et al.l2007l) and TRACE (jMcIntosh fc De Pontieu 
2009h time series. These upflows have been pro- 
posed as a source of the slow solar wind . Linear 



force-free (|Harra et al.l 120081 : IBaker et al.l 12 009) 



and potential field source surface (jSakao et al. 



tual uncertainties are of the order of 2 km, 
immediately sees that the sign itself of the veloc- 
ities (flows) may be changed. This clearly shows 
the need of a very careful determination of the 
wavelength reference. 

In this paper we set out to further explore the 
link between Active Regions 10942 and 10943. In 
Sec. |5] we use EUV image and spectroscopic obser- 
vations to analyze the flows in AR 10943. A spe- 
cial attention has been paid to the determination 
of the velocities taking account of the global flow 
velocities depen dance on temperature in the so- 
lar a tmosphere ( Chae et al. 1998 : Peter fc Judge! 
1999h . To understand the magnetic connection 



2007) extrapolations show open field lines in agree- 



ment with this hypothesis. 

Harra et al.l ( 20081 ) have also noticed that 



between both regions, we compute magnetic field 
extrapolations in Sec. [3] Finally, the results are 
discussed in Sec. 2] and we conclude in Sec. [5] 
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Fig. 1.— STEREO B/SECCHI EUVI observations of AR 10942 (on the left) and AR 10943 (on the right) 
taken at 195 A on 2007 February 20 at 8:05. The Hinode/EIS FOV at 05:47-07:59 UT for AR 10943 and at 
11:16-11:37 UT for AR 10942 for the data analyzed in this paper are drawn as rectangles. 
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Table 1: List of observations on 2007 February 20: Field-of-view (FOV), heliocentric center of FOV, start 
times (T s tart) or reference times (T re f), end times (T en( j) 5 and exposure times (T oxp ), center wavelength (A). 
Start and end times correspond to the period during which raster observations are used. 



Instrument 


FOV (arcsec) 


Center (arcsec) 


Tstart Or Tref 




T 

± cxp 


A or filter 


STEREO/SECCHI EUVI 


Full disk 




06:35 UT 


06:35 UT 


2s 


195A 


SoHO/MDI 


Full disk 




08:05 UT 




300s 


6767.8A 


Hinodc/EIS 


128 x 512 


(93, -37) 


05:47 UT 


07:59 UT 


60s 


195A 


Hinode/EIS 


41 x 400 


(-464, -100) 


11:16 UT 


11:37 UT 


30 s 


195A 
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Fig. 2.— Full-disk SoHO/MDI magnetic field measurements on 2007 February 20 at T re} = 8 : 05 UT 
clipped to the area used for the magnetic field extrapolations. Both AR 10942 (left) and 10943 (right) are 
visible. The Hinode/EIS FOV (at 05:47-07:59 UT for AR 10943 and at 11:16-11:37 UT for AR 10942, see 
Sec. 12. 2[) are drawn as rectangles. 
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2. Observations of AR 10943 and AR 10942 



2.2. Spectroscopic observations 



2.1. Observation with STEREO B/SECCHI 
EUVI image 

On 2007 February 20, less than four months 
after its launch, the heliocentric separation an- 
gle between the STEREO B probe and the Earth 
was still negligible (0.1 deg). This means th at 



STEREO B/SECC HI EUVI images (IKaiser et al. 



2008; Howard et al.l [20081) have the same viewing 
angle than SoHO and Hinode and that we can 
use these images in combination with SoHO/MDI 
and Hinode/EIS. STEREO/SECCHI EUVI was 
in its normal mode, with full-disk observations at 
a cadence of lOmin in the EUV channel A195A. 
We selected the 8:05 UT observation corrected by 
EUVI_prep from the SolarSoft library, shown in 

Figru 

In this image, both active regions can be seen 
simultaneously On the eastern side, the EUVI im- 
age in 195 A displays the AR 10942 complex made 
of a set of loops connecting the two extreme polar- 
ities (see Fig. [2]) of the AR, mainly in the south- 
ern side. Some straight (mostly fan-like) struc- 
tures are also clearly seen on the north-eastern side 
(X = -400, Y = -50 to 0, see Fig.[TJ) which are can- 
didate for open magnetic fields. At South, below 
AR 10942, internal loops, rather compact struc- 
tures (X = -300, Y = -300) seem to be the feet of 
(apparently) very sheared loops. On the western 
side of the image, another smaller and compact 
AR (10943) does not seem to be connected with 
its neighbouring regions, except for a set of diffuse 
loops at the East of AR 10943 whose feet seem to 
be located in between the two ARs. One also notes 
that the lower side of these diffuse loops is very 
sharp. The overall picture is that the two ARs 
seem to be very disconnected on one hand, and 
that the eastern feet of the diffuse loops mentioned 
above co uld coincide with t he QSL labelled "e" in 
Fig. 4 of iBaker et~aH (|2009l ). on the other hand. 
The EUVI image in the hot 284A line confirms 
this picture, contrary to the He II image where 
the chromosphere between the two AR does not 
seem to be very perturbed. Finally, the 171 image 
is more puzzling because it does not display the 
above (too) hot loops but also shows some dark 
area on the western side of AR 10942, which could 
be a coronal hole or a filament channel. 



Spectroscopic information can be obtained from 
the Hinode/EUV imaging spectrometer (EIS; 
Culh ane et al. 2007). We select two raster scans 
on 2007 February 20 that covers part of AR 10943 
(at 05:47 UT, study ID 45) and part of AR 10942 
(at 11:16 UT, study ID 57) respectively in order to 
have full spectra (around selected spectral lines) 
as a function of both solar dimensions in this re- 
gion. The slit positions during these raster scans 
are shown in the STEREO image (Fig. [1} and in 
TableQ] Both raster scans are partial on the active 
regions. Nevertheless, as we show with magnetic 
field extrapolation in Sec [3] the feet of the inter- 
connecting loops between the two active regions 
are located at the respective edges of the active 
regions in the FOV of EIS. So the raster scans are 
sufficient for our study. The delay of the scans 
and between the scans are negligible in compari- 
son with the time of the continuous flows which 
are visible on a few days. 

We apply the correction procedures eis_prep 
and eis_slit_tilt from the SolarSoft library. 
An additional correction must be applied for the 
orbital temperature variation of EIS; for this we 
have chosen to develop a specific method which 
is described in Appendix [A] : we use the orbital 
variation in the He II 256.32 A line to correct the 
orbital lineshifts in other lines. Indeed, this line is 
chromosphcric and optically thick, the insensitiv- 
ity of its centroid with respect to activity allows 
us to better isolate orbital variations. 

We focus on the Fe XII 195.12 A line which is 
emitted around logT = 6.1. We produce intensity 
(Fig. EJi for AR 10943 and Fig. [§i for AR 10942) 
and Doppler velocity (Fig. [5]d for AR 10943 and 
Fig. H)d for AR 10942) maps deduced from the pa- 
rameters of a single Gaussian fit of this line using 
the correction of orbital variation. We choose to 
ignore the self-blending of Fe XII 195.12 A with 
Fe XII 195.18 A because we do not focus on the 
width (which is the most influenced parameter) 
and our study concerns low density structures 
where the contribution of Fe XII 195.18 A is neg- 
ligible (|Young et al.ll2009l) . 

The level 2 data show that the velocities in 
the active region are of the same order in the 
other EIS windows available except for Fe XIII 
196.54 A where that pattern is reversed. In the 
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Fig. 3.— Maps of AR 10943 observed on 2007 
February 20 at 05:47-07:59 UT with Hinode/EIS 
in Fe XII 195. 12 A: a) Intensity (integrated on 
the wavelength dimension); b) Centroid wave- 
length; c) density from the Fe XII 196.64 A/ Fe XII 
195.12 A line ratio; d) Mass flow through a surface 
perpendicular to the line of sight. 



Fig. 4.— Maps of AR 10942 observed on 2007 
February 20 at 11:16-11:37 UT with Hinode/EIS 
in Fe XII 195.12 A: a) Intensity (integrated on 
the wavelength dimension); b) Centroid wave- 
length; c) density from the Fe XII 196.64 A/ Fe XII 
195.12 A line ratio; d) Mass flow through a surface 
perpendicular to the line of sight. 
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Fe XII 196.64/195.12 i 



Fig. 5. — Density function of the Fe XII 

196.64 A/Fe XII 1 95.12 A intensity ratio accord- 
ing to CHIANTI (iDere et all 12009} ) . The dashed 
lines represent the average of the intensity ratio 
and the corresponding density. The dotted lines 
represent the error bars on the intensity ratio com- 
puted from the uncertainties on the measurements 
of intensities and the related density error bars. 

Fe XII 186.88 A and Ca XVII 192.82 A lines, the 
core of the AR is blue but the other structures 
are the same. We cannot conclude about the 
Fe XIII 196.54 A and Fe XII 186.88 A lines veloci- 
ties because they are far different from the Fe XIII 
202.04 A and Fe XII 195.12 A respectively. The 
automatic analysis could not be sufficient. The 
Ca XVII 192.82 A line shows that the flow is up- 
ward in the core of the AR at very high tempera- 
ture (log T max = 6.7) but there are still downfiows 
in the vicinity of the active region. 

However the flow rate through a surface perpen- 
dicular to the line-of-sight is higher there than in 
the rest of the whirl. These downfiows are persis- 
tent for a few days, therefore an impulsive event 
such as a jet is excluded as a viable mechanism 
producing these localised redshifts. 

One also notes a redshifted area in the top third 
of the raster FOV. This area is not studied here for 
two reasons: it is not magnetically connected to 
AR 10942 (see Sec. |OJ) and the redshift is rather 
weak (a few km.s -1 ). 

We get a density map (Fig. \& for AR 10943 
and Fig. SJ: for AR 10942) by also fitting the Fe XII 
196.64 A line, and computing the Fe XII 196.64 A/ 
Fe XII 195.12 A intensity ratio, which is sensitive 



to density. The density is deduced from the theo- 
retical intensity ratio produced by the CHIANTI 
atomic database ( Dere et al . 1997, ,2009), as shown 
in Fig. [5l 

In order to derive the flow rate through a sur- 
face perpendicular to the line-of-sight (Fig. [3jl for 
AR 10943 and Fig. EH for AR 10942), we multiply 
the density by the Doppler velocity. 

The core of AR 10943 is characterized by hot 
loops bright in intensity (with a maximum of 9500 
counts/pix, Fig. [3^), downward Doppler velocities 
and high densities between 4 x 10 9 cm~ 3 and 1 x 
10 10 cm -3 . These densities are consistent with val- 



ues in active region loops by IWarren et al.l (2008) 



(1.3 x 10 9 and 9 . 5 x 1 10 cm- J for Fe XII) and 



with lYoung et al.l (|2009h (3 x 10 8 cm" 3 < n e < 1 x 
10 11 cm -3 ). The whirl of faint plasma around the 
core of AR 10943 is mostly blueshifted, with elec- 
tron density between 6.3 x 10 s and 1 x 10 9 cm~ 3 ; 
moreover there is a clear straight redshifted (up 
to 16 krn.s" 1 ) structure cutting the whirl in the 
South-East edge where the density is notably low 
(around 5 x 10 8 cm~ 3 ) but higher than Quiet Sun 
densities for Fe XII (n e = 2.5 - 3.2 x 10 8 cm~ 3 , 
I Warren fe Brooksll2009l) . 



3. Magnetic connection between AR 10943 
and AR 10942 

3.1. Magnetic field observations 

We investigate the possible magnetic connec- 
tivity between the two active regions AR 10942 
and AR 10943. We use a SOHO/MDI level 1.8 
96 minutes line-of-sight magnetogram (see Fig(2} 
to study the distribution of the photospheric mag- 
netic field. The SOHO/MDI magnetograms have 
been recorded on 2007 February 20 at T re f = 
08 : 05 UT (between the raster scans analysed in 
Sec. H] and simultaneous with the STEREO SEC- 
CHI/EUVI image of fig. [T]). We select an area en- 
compassing the active regions (heliocentric coordi- 
nate x from -720 to 275 arcsec and y from -321 to 
177arcsec). The total unsigned magnetic flux for 
this area is 4.05 x 10 22 Mx and the flux unbalance 
is only about 1.7% (a very low value for extrap- 
olation). The total unsigned flux for AR 10942 is 
9.90 x 10 21 Mx with a negative flux of 4.44 x 10 21 
Mx and a positive flux of 5.46 x 10 21 Mx. For 
AR 10943, the total unsigned flux is 5.53 x 10 21 
Mx with a negative flux of 4.12 x 10 21 Mx and 
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a positive flux of 1.41 x 10 21 Mx. The net mag- 
netic flux for AR 10943 is about 50% of the total 
flux. AR 10943 is in excess of negative flux, while 
AR 10942 is in excess of positive flux. If a mag- 
netic connection exists between the two active re- 
gions, thus this connection is between the positive 
polarity of AR 10942 and the negative polarity of 
AR 10943. 

3.2. Magnetic field extrapolation 

To study the connectivity between the two ac- 
tive regions, we derive the 3D magnetic fi eld based 



on a potential field extrapo lation (e.g., ISchmidt 



1964ISemel k Ravrole|[l968l ). The potential field 
assumption cannot describe the effects of elec- 
tric current on the magnetic field lines such as 
shear arcades or twisted flux bundles. However, 
the potential field assumption is a good estima- 
tion of the large-sc ale connectivity of t h e mag- 



netic field (see.g., Schatten et al. 19691; Levine 



1982 



Wang fc Sheelevlll992t ). In addition. |S gmer 



(120111 ) showed that, by comparing several mag- 
netic field assumptions (potential, linear force- 
free, and nonlinear force-free fields), the magnetic 
topology and connectivity are well preserved when 
the current density is modified, which justifies the 
use of the potential model in this study. 

The potential field extrapolation method solves 
the following equation: 



VxB = 0, 



(1) 



with the bottom (or photospheric) boundary con- 
dition given by the vertical component of the mag- 
netic field provided by the SOHO/MDI magne- 
togram selected area (see Sec. l3.1|) . The numerical 
technique relies on the computation of the scalar 
potential associated to the magnetic field satisfy- 
ing the Laplace equation. We perform the poten- 
tial field extrapolation for closed boundary condi- 
tions for which the normal component of the mag- 
netic field vanishes on the sides and top bound- 
aries. For the use of the SOHO/MDI magne- 
togram as boundary conditions, we need to put 
constraints: 

(i) the center of the FOV is near the disk cen- 
ter and the active regions are in a range of 
longitudes between E30 and W10, therefore 
the line-of-sight magnetic field is not subject 
to projection effects; 



(a) 



-400 -200 

Solor X (orcsec) 



Closed, z=10Mrr 



(b) 




(c) 



Fig. 6. — Connectivity maps for the potential field 
model in the case of closed boundary conditions, at 
two different altitudes: at the photosphere z = 
(a) and at z = lOMm (b). The color in any given 
pixel represents the length of the closed magnetic 
field line starting from this point (at the given al- 
titude) and the modulus of the magnetic field at 
this position, according to the caption in panel (c). 
The magnetic field is color-coded from 10 to 100 G 
and saturated below and above these values, while 
the the length is coded from 10 to 300 Mm and sat- 
urated below and above these values. Rectangles 
in (b) indicate the EIS FOV. 
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(ii) the line-of-sight magnetic field component 
B s is converted into the vertical (radial) 
component B z — B s j cos 9 where 9 is the 
angle between the line-of-sight and the nor- 
mal to the surface. 

(iii) as the magnetic field component satisfies the 
two items above, we assume that the cur- 
vature of the Sun does not affect the con- 
nectivity of the field. We thus compute the 
potential field in Cartesian coordinates. 

(iv) as the magnetic flux through the surface con- 
taining the two active regions is nearly bal- 
anced (<2%), we do not force the magnetic 
flux to be strictly balanced. 

In Fig. we plot the connectivity maps for the 
potential model with closed boundary conditions: 
each pixel is color-coded as a function of the length 
of the field line associated with the photospheric 
footpoint and MS c\ function of the magnitude of 
the field, according to the caption in panel (c). 
The connectivity maps have been computed at two 
different altitudes: at the photosphere z = (a) 
and at z = 10 Mm (b). The modulus of the mag- 
netic field is below 10 G in the Quiet Sun, mostly 
between 10 G and 100 G in the active regions at 
10 Mm while superior to 100 G in some active ar- 
eas at z = 0. Long field lines linking the two 
active regions definitely exist, when starting from 
both altitudes. At the photospheric level, we no- 
tice that the connectivities have various sizes and 
magnitudes, the smallest (and faintest) ones cor- 
responding to local dipoles in the quiet Sun. On 
the contrary, at z = 10 Mm, these local dipoles 
disappear, probably because of their small size, 
and consequently altitude. One recognizes the es- 
sential long range connectivity between the two 
AR, along with an internal connectivity within 
AR 10942. This large scale connectivity between 
the two A R has already been found using different 
models by Sakao et al. (120071) . who did not discuss 
it, and by lHarra et al.l (|2008l) . 

3.3. Estimation of the mass flows 

We explore further the magnetic field connec- 
tion between the ARs, using the potential field ex- 
trapolation (Fig.[7K and b) and EIS spectral scans. 
The co-alignemcnt between maps of EIS inten- 
sity and maps of the extrapolated magnetic field 



is made by using cross-correlation between MDI 
(8:05 UT) and STEREO Fe XII 195 (8:05 UT) im- 
age and cross-correlation between STEREO Fe XII 
195 image and EIS Fe XII 195.12 (respectively 
at 05:47 UT for AR 10943 and at 11:16 UT for 
AR 10942). 

For each AR, we have measurements of the line- 
of-sight velocity v s (oriented downwards) , and the 
signed mass flux through a surface S orthogonal 
to the line-of-sight is 



F = 



pv ■ d 2 S 



(2) 



which reduces to pv s S if S is small enough. Note 
that F has the same sign as v s . 

We choose to consider the two sections Si and 
S% of the flux tube represented in Fig. [5] We de- 
compose S\ into small Su and then we define 
Sai by the intersection with 5*2 of the flux tube 
originating from Su: we have Si = Y). Su and 

S2 = Y^i S2i- 

The ratio of the total mass flux through S2 to 
the total mass flux through Si is 



r = Ei = ^iEii 



Yi PliVsliSn 
Yi Pli v sli^li 



(3) 



We use the conservation of the magnetic flux 1 
in each of these small flux tube^|: 



10* 



\B s ii\Su — \Bg2i\S2i 



(4) 



and we decompose into SuS^i/ 'Su, which gives 
S21 = Su\B s u/ B S 2i\- Finally, as all Su are chosen 
with the same area (a MDI pixel), we obtain the 
mass flux ratio 



R 



Yli P2iV s 2i\B s u/B s 2i 
Yli PliVsli 



(5) 



In this equation, the density and the Doppler 
velocity from the EIS Fe XII spectra in AR 10942 
and 10943 (see Sec. 12. 2|) are derived with a 3 x 3- 
pixel median filter. We consider that the Fe XII 
lines are emitted at an altitude of 10 Mm in hot 
loops. The values of the line-of-sight magnetic 
fields Bu and B 2 i are taken from the footpoints 



x Flux and magnetic field are in absolute values because the 
surfaces' orientation is given by the line-of-sight instead of 
the direction of B. 
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Fig. 7. — Some selected magnetic field lines from a potential magnetic field extrapolation, with the MDI image 
as background: (a) Viewed from above with superimposition of the velocity maps from EIS for AR 10942 
and AR 10943; (b) Viewed from a point of view in the Southern hemisphere; (c) Start (black crosses) and 
end (white crosses) points of the field lines used in Sec. !3.3l on a map of B z at the 10 Mm altitude obtained 
from the extrapolation; three selected lines are drawn in black. The black frames arc the EIS FOV. 
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of magnetic field lines starting at an altitude of 
10 Mm on a 7 x 15-point grid in the outflow 
at the edge of AR 10942 (Fig. Efc), and com- 
puted until they reach again the same altitude of 
10 Mm. 76 among 105 of these lines link AR 10942 
and AR 10943 and fall inside the EIS FOV on 
AR 10943. 

We consider these lines to be representative of 
the flux tubes connecting AR 10942 to AR 10943. 
The modulus of the mass flux ratio (AR 10943 over 
10942, see Eq. [3) is then 0.18. 



and 




Fig. 8. — Cartoon of the magnetic flux tube 
connecting AR 10942 and AR 10943. Si and S 2 
are the cross-sections perpendicular to the line- 
of-sight of the flux tube at the observed altitude 
(we chose 10 Mm for Fe XII), above AR 10942 
and AR 10943 respectively. The magnetic field in 



cross-section Si is Bi, 
velocity is Vi. 



the density is pi, and the 



We compute the uncertainty on the mass flux 
ratio induced by the uncertainties on the observa- 
tions by 

SR _ 5F 1 
R ~ \FA 



6F 2 



(6) 



with 

5 Fl _£E,l*i<l{ft 



1*11 



\Fi\ 



1 I I ( Sfti , $ v sli 

-Pi K/n ^ \ pu \v s ii\ 



(7) 



SF 2 



^ Hi 1*2*1 \p£\ 



\Fi 



F 



7- 1 =y^P2i\v S 2iB sli /B s2i \x (8) 
6B s u 5B s2 i 



5p 2 



Sv x 



\Vs2i 



\B aU \ 



\B S 



We have chosen a uniform Sv of 2 km.s" as we 
have an uncertainty of about zblkm.s -1 for the 
Doppler velocity from each line fit, Sp — 3.4 x 
10 8 cm~ 3 from the uncertainties on the intensities 
of the lines used to deduce the density (FigEJ), and 
SB = 10 G, a value below which coherent small- 
scale patterns disappear in the quiet Sun. The 
relative uncertainty 5R/R on the mass flux ratio 
is then found to be 0.30. 

To evaluate the uncertainty due to the co- 
alignment between MDI data and EIS spectra, we 
compute the mass flux ratio with a co-alignment 
shifted by one MDI pixel (1.98arcsec) in the two 
X and Y directions with respect to EIS data. The 
mass ratios obtained vary by a factor 2. We can 
then state that the co-alignment is the main cause 
of uncertainty. 

4. Discussion 

The extrapolation clearly shows us how AR 10942 
and AR 10943 are linked. From the observation 
of both ARs, it appears that the mass flow to 
AR 10943 seems to be 18% of the upflows from 
AR 10942. We can therefore conclude that the flow 
leaving the Sun among the upflows from AR 10942 
is at most 82% of the measured upflows. Some of 
these upflows could be confined somewhere else. 

Another solution is not explored by our extrap- 
olation: the presence of a Y point between the 
two ARs. An argument in favour of this hypoth- 
esis is the difference between the velocities at the 
AR 10942 edge and the AR 10943 edge. I ndeed , 
in the model proposed by iDel Zanna et al.l (120111 ) , 
the velocity is higher for the strongest magnetic 
bipoles than in the faintest pole (a Coronal Hole 
in their model). In this case, the flows are not 
directly linked and they interact with each other 
through an X point. In our case we clearly have 
open field lines, with the presence of two coronal 
holes, one located at the Eastern part of AR 10942 



11 



and — with less evidence — the other one in the 
Western part of AR 10943, which then would lead 
to a high-altitude Y point (Fig. [9|). As a conse- 
quence, the measured upflows in the Eastern part 
include flows of material which escapes towards 
the Y point and the corona and cannot be found 
at the foot of the connecting loop in AR 10943. It 
is clear that the potential extrapolation does not 
allow to explore this possibility. 



5. Conclusion 

In summary, we observed Doppler-shifted struc- 
tures at the edges of two Active Regions which 
we find to be connected with each other by large 
scale loops, a result we discuss in more details be- 
low. With the advent of the EUV spectrometer 
EIS on Hinode, special attention has been paid 
on high temperature flows around coronal holes 
and active reg ions. AR 10942 ha s been extensively 
studied (e.g. Baker et al. 20091) and it has been 
suggested that the strong observed upflow at its 
boundary could be one of the sources of the solar 
wind. In order to properly assess this critical re- 
sult, we have investigated the possibility that the 
supposedly open field in this area could be actu- 
ally (at least partially) connected to the preced- 
ing AR 10943 with the consequence that part of 
the mass outflow from AR 10942 could turn into 
downflow in AR 10943. In order to investigate this 
possibility, we analyzed the EIS spectra obtained 
in AR 10943 with special attention to the Fe XII 
(195 A) line formed at logT = 6.1. 

The EIS spectra are difficult to analyze for two 
reasons: most lines are blended and their profiles 
are affected by a strong (presumably thermal- 
induced) deformation during any orbit. The 
blends are properly dealt with the use of multi- 
gaussian fits. As far as the second (orbital) correc- 
tion is concerned, the EIS software allows to take 
into account the relative shift of the profile over 
the orbit but leaves open the issue of the absolute 
reference. In order to overcome this difficulty, we 
have designed and used a technique which relies 
upon the cool and optically thick 256 A He II line, 
consequently less sensitive to activity and induced 
flows. The technique is described in Appendix lAl 

The analysis of the velocity field in AR 10943, 
performed with this technique, shows a significant 
feature located in the southwestern part where the 




CH AR 10942 



AR 10943 



Fig. 9. — Cartoon of the magnetic field and flows 
in the case of a Y point. The colors of the lines 
show the Doppler shifts that would be observed 
from above, with the usual convention: blue for 
upflows and red for downflows. The Y point is 
located above the colored loop and is shown with 
the current sheet in grey. Some field lines are in- 
dicated with arrows. At the bottom, the magnetic 
polarities coded in the usual white and black con- 
vention correspond to the location of the two ARs 
and the eastern coronal hole. 
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downflow velocity reaches 16km.s _1 . It should be 
pointed out again that this pattern is confirmed in 
other lines, as mentioned in Sect. 12.21 with the no- 
ticeable exception of the Fe XIII line at 196.5 A 
which is a very weak line ( Young et al. 2009). 



At very high temperature (Ca XVII 192.82 A), 
it seems that the pattern is the same except for 
the core of the active region which is blueshifted. 
Moreover, through the use of the two Fe XII lines 
(196.6 A and 195.1 A) and their intensity ratio, we 
were able to build a density map where the down- 
flow region density reaches a value of 5.10 s cm -3 . 

We then proceeded to establish the connectiv- 
ity between the two regions, and especially be- 
tween the blueshifted area of AR 10942 and the 
redshifted area of AR 10943. A potential extrapo- 
lation, rather suitable for our stable active regions, 
was used. The connectivity between the two ar- 
eas was well established although the computed 
field lines do not clearly appear as inter-connecting 
loops in the Fe XII image (Fig. [I}, presumably a 
result of their low densities. With the values of the 
local density and line-of-sight velocity in the two 
areas, we could compute and compare the mass 
flows between the respective feet of the magnetic 
field lines. We concluded that only 18% of the out- 
going flow from AR 10942 is retrieved as downward 
flow in AR 10943, as computed along the extrap- 
olated field lines. We elaborated upon the pos- 
sibility of other paths followed by the upflows of 
AR 10942, for various reasons: one of them being 
the presence of open fields in AR 10942, which our 
extrapolation does not take into account; another 
one being the possibility of a breach of connec- 
tivity d ue, e.g. to a Y point where reconnection 
occurs ( Del Zanna et al.ll2011 ). 

As a conclusion, we can state that some mate- 
rial is exchanged between the two active regions, 
especially from the area of the upflow but that 82% 
of the upflow from AR 10942 is missing upon ar- 
rival on AR 10943. Actually, these flows could be 
a condition for the thermal sustainability of such 
long loops during such a long time span. How- 
ever, we cannot go further because of the many 
factors and uncertainties involved in the observa- 
tions, whether the wavelength (velocity) reference, 
the non-simultaneity of the two AR observations 
with EIS, the reduced fields of view, or even the 
validity of the potential field extrapolation. 

We strongly recommend that further simi- 



lar studies rely upon force-free magnetic field 
extrapolation made possible by simultaneous 
SOT/Hinode observations and also closer-in-time 
observing modes focused on the active regions 
involved. We also suggest to use the reference 
wavelength method described in the Annex. 
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A. Wavelength reference for EIS 



The temperature varia tion along the orbit of Hinode deforms the optical axis of EIS in the dispersion 
plane (jBrown et al.ll2007|) . The determination of the Doppler velocity is therefore more difficult because of 
the lack of absolute reference. Figure fTUl shows the wide range (0.05 A) of the Fe XII wavelengths giv en by 
EIS before any correction. Therefore a correction is necessary. It has been observed (jBrown et al.ll2007|) that 
this orbital variation of the wavelength is not the same from one orbit to the other and depends on several 
parameters (Hinode pointing, Earth's seasons, etc.). Thus, it is not predictable and the orbital correction, 
based on the considered scie nce data set, has to be empirical, although a correction where housekeeping data 
are used ( Kamio et al.ll20foT) is possible. Nevertheless, as the orbital variation is t he same for both detectors, 
we can use the observation of any line to correct the others (jMariska et al.l 120071 ). 

The measured wavelength depends on the physical parameters on the Sun and the instrumental bias 
including the orbital variation. We want to separate these two components of the wavelength line centroid. 
We assume that the Doppler shifts are globally null on average in the Y dimension while the orbital variation 
is the same on each pixel along the Y dimension of the slit at a given time. So we can isolate the orbital 
variation by integrating on the Y axis. In order to smooth any residual Doppler shifts, we use a spline 
depending on time. 

The procedure consists in averaging the line profiles along the Y axis in the quiet Sun pixels and then 
to fit these integrated lines with a Gaussian in order to obtain the mean wavelength of the line centroid 
depending on the time X av (t). The line centroids X av (t) obtained are fitted with a spline S(t)) in order to 
keep only the orbital variation. The correction, that is to say the difference between the spline and the 
reference wavelength Ao (as given in CHIANTI), is denoted 5X(t). 

The usual way of finding the wavelength reference consists in following the above procedure. However, the 
issue of the line shift in the corona is not solved today. From SUMER/SOHO, it was derived that the redshift 
observed in the transition region (with a maximum of 30 k m.s" 1 at logT= 5.3) decreased at high coronal 
tempe ratures but remained a redshift (e.g. Chae et aL"lll998 ). Then also with SUMER data, Peter fe Judgel 
(1999) showed that at high temperatures, one actually had a blueshift, i.e. upward velocities. Consequently, 
we suspect that taking a hot line as a wavelength reference, even with a substantial spatial averaging, may 
not be a safe procedure. 

As the He II 256.32 A line is the coldest line in the EIS wavelength ranges (being formed from the upper 
chromosphere to the lower transition region) and is optically thick, the line is less affected by the solar flows 
than the other EIS lines. Therefore, the centroid of its profile is an efficient reference. The drawback is 
that He II 256.32 A is blended wi th hot lines. In the quiet Sun, He II 256.32 A is dominant (its contribution 
is more than 80% of the blend (jYoung et al.l 120071 ) and can be separated from the blending lines with a 
two-Gaussian fit. For the active Sun, it is more difficult to separate He II 256.32 A from its blends, as they 
are hot lines and their contribution is changing depending on the observed structures. 

In our study, for example, we choose to apply the correction obtained from He II 256.32 A on the Fe XII 
196.64 A line. We computed 5\Heii(t), we fitted \FeXii(y,t) by a gaussian fit and we corrected each pixel 

following the formula X Fe XIIcorrected(y,t) = X FeX n(y,t) - 8X H eii{t) 

As He II 256.32 A is essentially blended with Si X 256.37 A (|Del Zanna fe Masonl liooBh . We consider 
that we can fit these two lines only. Because Si X 261.04 A and Si X 256.37 A have the same lower level 
(2s2p 22 Pi/ 2 ), we can use the branching ratio between Si X 261.04 A and Si X 256.37 A in order to evaluate 
the contribution of Si X 256.37 A in the 256 A blend. The branching ratio denoted a is determined by the 
ratio of the intensities given by CHIANTI. 

The procedure is as follows : 



• At first, we fit Si X 261.04 A with a simple Gaussian. We obtain XsiX26i, 1 Si X26i, FW H M si X26i an d 
the background. 



We transpose the characteristics we obtained for Si X 261.04 A to Si X 256.37 A, that is to say for 
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the wavelength. X Si X 256 = ^SiX26i ~ Ao,siX26i + ^o,SiX2m, we consider the same width, and for the 
intensity, I256 = &I261 where a is the branching ratio. 

• We fit He II 256.32 A and its blend with a 5-parameter fit (\hcII- IhcII- Isi X256> FWHMu e ii an d 
the intensity background), with fixed value for \si X25& and FWHMsiX256- 

• We extract the correction SX(t) from Ajj e j/(t) as described previously and subtract it from the wave- 
length obtained for the studied line. 

Figure [HI shows the dispersion of the wavelength of Fe XII 195.12 A after the correction proposed in this 
appendix. We can note that the wavelength increases slightly with the intensity, a result in agreement with 
the fact that the Active Region is slightly redshifted. The orbital correction computed with the SolarSoft 
library based on the Fe XII 195.12 A line shows the same tendency (Fig|l2]). There is a systematic difference 
of 0.008 A b etween the t wo co rrectio ns which is cons i stent with the line shift depending on the temperature 
discussed in lChae et al.l (1l998 ) and i n Peter fe Judge ( 19991 ). They reported a redshift of 5.3 km.s -1 for Si III 
A 1206.51+A 1206.53 A (jChae et al. l ll998l) at ab out the same temper ature as He II 256 (\ogT max = 4.7) and 
a blueshift of 6 km.s" 1 for Fe XII 1349 A line l)Sandlin et al.lll977t) which corresponds to a total relative 
Doppler shift of 0.008 A. Therefore, we apply a redshift of 0.008 A on our Doppler shifts analysis (Sec 
and Sec. 13. 3p . which implies that we adopt an average quiet Sun Fe XII line line at rest. 
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Fig. 10. — Dispersion of the centroid wavelength of the Fe XII line profile depending on the intensity (counts) 
before the orbital correction. 
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Fig. 11. — Dispersion of the centroid wavelength of the Fe XII line profile depending on the intensity (counts) 
after the orbital correction using He II. 
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Fig. 12. — Dispersion of the centroid wavelength of the Fe XII line profile depending on the intensity (counts) 
after the orbital correction using Fe XII. 
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